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1. Introduction 

The TAUOLA package jll2l3l4j for the simulation 
of r-lepton decays and PHOTOS |5I6| for the simu- 
lation of radiative corrections in decays, are com- 
puting projects with a rather long history. Writ- 
ten and maintained by well-defined authors, they 
nonetheless migrated into a wide range of appli- 
cations where they became ingredients of com- 
plicated simulation chains. As a consequence, a 
large number of different versions are presently 
in use. From the algorithmic point of view, they 
often differ only in a few small details, but in- 
corporate many specific results from distinct r- 
lepton measurements. Such versions were mainly 
maintained (and will remain so) by the experi- 
ments taking precision data on r leptons. On 
the other hand, many new applications were de- 
veloped recently, often requiring a program inter- 
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face to other packages (e.g. generating events for 
LHC, LC, Belle or BaBar physics processes). 

In the last two years some progress in the sim- 
ulation of r-pair production alone was achieved. 
A control of the systematic errors was improved 
and a special tool developed for that purpose, 
MC-TESTER, turned out to be very useful. With its 
help, the new version of TAUOLA, which may be- 
come the starting point for the work directly ori- 
ented towards phenomenology of b-factories, was 
prepared. MC-TESTER became instrumental in de- 
veloping the new version of PHOTOS for radiative 
corrections in decays. With its help tests enabling 
verification of the program in multiple photon ra- 
diation mode became possible, thus establishing 
PHOTOS as a precision tool for simulation of ra- 
diative bremsstrahlung in W, Z and Higgs boson 
decay. New results for r physics were obtained as 
well. 

Because of the limited space and sizeable 
amount of physically interesting numerical re- 
sults, only some of them will be included in the 
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conference proceedings. For the remaining ones 
we refer the reader to references and transparen- 
cies from the talk. 

2. Versions of TAUOLA Monte Carlo 

In refs. |4I7| the setup for constructing specific 
versions of TAUOLA and PHOTOS from the single set 
of files was prepared and documented. The sys- 
tem was prepared for the software librarians and 
advanced users interested in updating the pack- 
ages for the multipurpose environment. The idea 
was to create a repository which allows one to 
keep all main options of TAUOLA developed for 
different purposes in a relatively compact form, 
without duplications of semi-identical parts. The 
repository was set to produce standard FORTRAN 
files which can later be handled exactly the same 
way as the early 1990 versions of the packages, 
but with physics initialization based on Cleo and 
Alcph data of the years 1997 - 98. During the 
conference, another generating system was de- 
veloped in parallel, to start new work with the 
Belle and BaBar collaborations. Physics initial- 
ization remained as in Cleo, but 60 new channels 
were prepared (without physical initialization) for 
the future implementations to be done within 
those collaborations. To facilitate integration, 
the generic subdirectory TAUOLA/tauola-BBB was 
prepared. The new versions can be now initial- 
ized with the help of the command make with 
parameters authors. Belle or BaBar. Whether 
this solution will become accepted needs to be in- 
vestigated. The geographical spread of interested 
people over three continents is certainly a disad- 
vantage. 

3. PHOTOS 

Recently some extensive tests and extensions 
were introduced into PHOTOS, an algorithm for ra- 
diative corrections in decays of particles or reso- 
nances. Thanks to comparisons with the single- 
photon exact matrix-element for the W decay cor- 
recting weight was introduced (in this particular 
case, see fig^l. The method applied can also serve 
as an example for similar correcting weights to be 
introduced in other decay channels, of more im- 



portance for the b-factories community. In theses 
cases correcting weights may be of a more phe- 
nomcnological nature, and to some degree may 
represent a fit to the data. 

In the second part of the extension of PHOTOS 
we have followed the solution presented in ref. 
[H] for the iterative algorithm of multiple photon 
emission. Thanks to a better understanding of 
the numerical stability of the algorithm, we were 
able to go beyond double photon-emission, and 
extend the iteration up to quadruple emission. 
Once this was achieved we were able to complete 
the exponentiation programme, still preserving 
full control of terms providing leading logs as well 
as the soft photon region of the photon emission. 
It is important to note that the algorithm covers 
full phase space for photon emission, except soft 
regions where photon emission is integrated an- 
alytically and added together to the virtual cor- 
rections. Thanks to the exponentiation, we could 
lower the value of the minimal energy of gener- 
ated photons to the level of 10~^ or 10~^ without 
any problems. This was plaguing the project un- 
til now. 

This goal required an important development 
for numerical stability. PHOTOS operates on four- 
vectors, that is why explicit generation of pho- 
tons at small energies immediately leads to loss 
of the numerical stability of additional 6 orders. 
Therefore a special algorithm to prevent accu- 
mulation of rounding errors was developed. On 
the technical side let us remark that this spe- 
cial correcting routine, PHCORK, had to be called 
after every emission of individual photons (rou- 
tine PHOMAK. PHOMAK has to assure that energy- 
momentum is conserved in the decay branch and 
that four-momenta of all particles are precisely 
on mass-shell. 

Once the list of necessary modifications was 
completed we performed extensive tests using 
the method presented in the next section. This 
helped to estabhsh the precision of PHOTOS in the 
case of Z decay at the permille level. Comparison 
tests with KORALZ and KKMC support such conjec- 
ture. In the case of W decay it was confirmed by 
comparison with WINHAC [Tnj . 
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4. MC-TESTER and its extensions 
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Figure 1. Comparisons (ratios) of the SANC and 
PHOTOS predictions for the W decay, distribution of fi^v 
acolinearity angle, cos{Oacoi) is given. Ratio of the dis- 
tributions from the two programs is plotted. The dom- 
inant contribution to the distribution is from the non- 
infrared, non-leading-log part of the phase-space, i.e. 
the region which is most sensitive to non-leading con- 
tributions 

Upper two plots: Distribution with correction 
weight, plot taken from ref. 8 . Results with approx- 
imation installed in SANC matrix element provide addi- 
tional technical test. 

Lower plot: Distribution without correction weight 
plot taken from ref. [5| 



MC-TESTER is described in detail in ref. [TT] : 

the principle of this testing package is to analyse 
series of events originating from a Monte Carlo 
generator (or a system of several Monte Carlos 
combined together). In fact only a branch origi- 
nating from the decay of some intermediate res- 
onance or particles is examined. In a fully auto- 
mated way (thereby limiting manpower and also 
the risk of accidental errors) the output from the 
run is created. If a similar run for a different 
generator (possibly even written in different pro- 
graming language) is performed, then the two files 
with information collected by MC-TESTER can be 
compared and the result of the test can be pre- 
sented in a visualised form as a booklet of plots 
and tables with general information on the run, a 
list of the decay channels, their branching ratios, 
and SDP^ parameter. 

Such a test was very useful in the case of val- 
idation of TAUDLA. However, in cases where an 
arbitrary number of final-state particles could 
be present (as in the case of final-state QED 
bremsstrahlung) the test required an extension. 

There may be quite big differences in the way 
bremsstrahlung is generated by the programs for 
decays. In our test, we aim to develop a tech- 
nique where comparisons make physical sense and 
are automatic. To this end, we have to free our 
test from dependence on the technical parameters 
of the generators, such as a minimal photon en- 
ergy threshold, hence ambiguities due to the num- 
ber of soft photons generated, explicit choices for 
Lorentz frames of generation, etc. 

Independently from how events are generated, 
we define zero-, one- or two-photon topologies. 
We will call the event to be of 'zero photons', if 
there is no photon at all of energies larger than 
the Etest parameter of the test. The 'one-photon' 
event will have to have one (and only one) photon 
of energy larger than Etest ■ If there is more than 
one such photon, we shall call it 'two-photon' 



^ Shape DifTerence Parameter: quantifies how much the 
shape of the histograms (of all possible invariant masses 
which can be constructed from the momenta of the decay 
products of the particle under study) differs between the 
two runs of distinct programs. 
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event. If there are more than two photons of 
energy larger than Et^st, then we consider only 
the two most energetic photons, and treat the re- 
maining ones as if they did not pass the Etest 
threshold. 

There may be cases where more than two pho- 
tons actually considered (e.g. with energies not 
passing the Etest threshold) are present. Fol- 
lowing the leading-log-inspired logic, we will add 
them to the momenta of outgoing fermions (al- 
ways to the fermion of a smaller angular separa- 
tion). 

We define two tests: testl and test2. Test2 
is exactly as explained above. In case of testl, 
only one photon (the most energetic one) can be 
accepted. 

The test has therefore two parameters: 

- maximal multiplicity (of photons) defines 
whether testl or test2 is used 

- minimum photon energy: Etest 

Tests of PHOTOS in W and Z decays 

Tests of comparisons for PHOTOS with KKMC and 
KORALZ |12ll3j for Z decay and with WINHAC HO] 
for W are collected on the webpage ^3]. These 
tests are quite extensive and we can not present 
here even the most important results. That is 
why only an example illustrating the general form 
of the results is given in fig. El The upper part of 
the figure presents the comparison 'hyper-table'. 
Comparisons of PHOTOS running with different 
options (single- double- triple- quadruple- and 
multiple-radiation with the results from KKMC and 
KORALZ Monte Carlo programs are collected. The 
first column of the table defines the generator 
used (with options) for the numerical results col- 
lected; the second part of the table collect branch- 
ing ratios obtained from testl and test2 as de- 
scribed above, respectively for configuration with 
0,1 (testl) and 0,1,2 photons (test2). Later the 
maximum of the SDP (as defined in ref. ^J) is 
given. As reference point, the first entry after 
the definition of Etest is used; in this table re- 
sults from KKMC are used. The last two columns 
represent hyper-links to the comparison booklets 
(not included in the articles but available from 
[H]). The booklets of standard MC-TESTER for- 



mat TT" include the plots like the one presented 
in the middle of fig.[21for invariant mass of ^jL~ 
where the largest difference between results from 
the KKMC and PHOTOS version with exponentiation 
is visible. 

One can observe, rather unexpectedly, a good 
agreement between PHOTOS and KKMC. In fact only 
in the case when the full double-photon matrix- 
element in KKMC is used is the agreement with 
PHOTOS good. If a less powerful matrix-element 
that does not rely on the spin amplitude expo- 
nentiation scheme of KORALZ is used, or a second- 
order matrix-element is switched off in KKMC, dif- 
ferences from the full exponentiated version of 
PHOTOS, quantified with the help of SDP, increase. 
What could be even more unexpected is that the 
effects of coherence between photons is to a large 
degree reproduced by the PHOTOS algorithm. This 
delicate point is still under investigation. Some 
numerical results are collected on the web page 
[H] in subsection 'PHOTOS improved antenna 
algorithm'. However, this point is of no numeri- 
cal interest for bremsstrahlung in decays, at least 
not for present-day applications, because the ob- 
served effects are too small. On the other hand, 
this point is interesting as a potentially important 
ingredient for constructing high-precision Monte 
Carlos for future applications; either for QED at 
linear colliders, or for QCD at the LHC. 

In the case of W decays our comparisons for 
multiple photon emission were limited to com- 
parisons with WINHAC lOj. Less good agreement, 
but still better than required by present-day ap- 
plications, was achieved. We did not investigate 
the sources of the residual discrepancies. The 
good candidates were limitations of the correcting 
weight used in PHOTOS, see discussion above, and 
the fact that the WINHAC matrix-element is limited 
to the first order only. More detailed numerical 
discussion of this point could not be included in 
the conference contribution presented because of 
space limitation. 

Tests of TAUOLA and PHOTOS in leptonic t 
decays 

In the case of the t lepton physics, the discus- 
sion presented above may be of no practical con- 
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sequences. However, the methodology developed 
for QED bremsstrahlung (and tests) in case of W 
and Z decays could also be applied for leptonic 
r decays. With the methods analogous as pre- 
sented above, we have found that the difference 
between single-photon full matrix- element gener- 
ation of TAUOLA and the single-photon option of 
PHOTOS is larger than between the first-order and 
exponentiated mode of operation of PHOTOS, see 
ref. We could conclude that the matrix el- 

ement of TAUOLA for r vv^{'-^) is better than 
the exponentiated version of PHOTOS. This conjec- 
ture is true at least for semi-inclusive quantities 
as discussed here. Of course for generation of 
events with explicit multiple-photon final-states, 
the results generated with PHOTOS will be more 
precise. 

5. Universal interface based on HEPEVT 
common block 

The universal interface of TAUOLA for 'any' r 
production generator is distributed together with 
TAUOLA [CTt] , It uses as an input the HEPEVT 
common block and operates on its content only. 
As a demonstration example the interface is com- 
bined with the JETSET generator [T^. 

The interface works in the following way: 

• The T-lepton should be forced to be stable 
in the r production generator. 

• The content of the HEPEVT common block 
is searched for all r leptons and neutrinos 
first. 

• We check if there are r-flavour pairs (two 
r-leptons or r-lepton and r-neutrino) orig- 
inating from the same mother. 

• The decays of the r-flavour pairs are per- 
formed with TAUOLA. Longitudinal spin ef- 
fects are generated in the case of the t pro- 
duced from the decay of: W tv, Z/7 — > 
rr, the neutral Higgs boson H tt, and 
the charged Higgs boson tv. Paral- 
lel or anti-parallel spin configurations are 
generated, before calling on the r decay, 
and then the decays of 100% polarized r's 
are executed. 



• In the case of the Higgs boson (for the spin- 
correlations to be generated) the identifier 
of the t's mother must be that of the Higgs. 
In this case it was rather easy to implement 
full spin-correlations. 

• In case the r lepton is not produced from 
any of the above-mentioned intermediate 
states, then if there is produced alsoa Vr 
from the same mother as that of the r , the 
W is reconstructed by the interface as the 
sum of the two. Similarly, Z/7 is recon- 
structed if another r is produced from the 
same mother. 

This organization of software became quite pro- 
ductive. It was helpful to perform studies of ob- 
servables for Higgs boson parity measurement at 
future LCs. In this case it was important that 
full spin correlations for the multi-layer cascade 
process h t~^t~, — > i^p*, — > tt'^tt^ be 
included |17ll8ll9l2fl] . Also, in the case of studies 
for the Higgs discovery potential at LC, the im- 
plementation of the universal interface was use- 
ful and productive to better visualize possi- 
ble difficulties. Recently [55], the interface had 
proven to be useful for studies of CP parity ef- 
fects in BqBo system produced at Belle/BaBar 
energies, when one of the B-mesons decays into 
a pair of r leptons. This process, even though 
not discovered so far, provides another example 
of the usefulness of the universal interface. 

6. Summary and future possibilities 

The status of the computer programs for the 
decay of r leptons and associated projects was re- 
viewed. The high-precision version of PHOTOS for 
radiative corrections was presented. In particu- 
lar, the option to run the program with multiple- 
photon radiation was presented. The results 
for decays of r-leptons, Z and W bosons were 
mentioned. It was found that for semi-inclusive 
single- and double-photon final-states the agree- 
ment with the matrix element simulations of the 
second-order is at the 0.1% level. Only at the tails 
of some distributions, the approximation used in 
PHOTOS was visible and gave a discrepancy of the 
15% level in sparsely populated corners of the 
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phase-space. 

The presentation of the TAUOLA general- 
purpose interface was ommited because of lack 
of time. Examples for its use in the case of the 
Higgs boson parity measurement at a future Lin- 
ear Collider |18I19I20| and for Higgs searches at 
the LHC can be found in the literature. Recently, 
a similar application was developed for the case of 
studies in hypothetical effects of CP-parity break- 
ing in the Bq-Bo system at Belle and BaBar |22|. 
At present, in its full content, it is available from 
the authors upon individual request only. 

Distinct versions of the TAUOLA library for t 
lepton decay, and of PHOTOS for radiative correc- 
tions in decays, are now in use. The principles of 
the distribution package, as presented in ref. |7j, 
are carefully preserved to save continuity of the 
project, and for the comfort of the users. In this 
scheme the new version of PHOTOS is now avail- 
able from the home page of Z. Was ^Sj. How- 
ever, on top of the previously available versions of 
TAUOLA, an extension was prepared to introduce 
new optional and/or missing decay modes from 
the independent files. In these generating files 
60 new channels of r decays were drafted on top 
of the previously available ones. As previously, 
the system for creating the required version of 
the PHOTOS and TAUOLA packages from the single 
master copies is preserved. The master copies are 
kept in relatively compact and clear form, with- 
out unnecessary code duplications. 

Development of MC-TESTER was a necessary 
step to prepare the translation of TAUOLA to new 
programming language such as C++, however, it 
is also very important to perform such transfor- 
mation together with the program's users. At 
present, there is no pressure from the experimen- 
tal community for such transformation to be im- 
minent. As far as the package authors are con- 
cerned, such a possibility has been envisaged for 
some time now. 
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Figure 2. Upper part: Comparison table for 
PHOTOS, KKMC and KORALZ results. For more details see 
the text. Middle plot: Single plot from the compari- 
son booklet. The chosen plot presents the distribution of 
the ^^/i^ invariant mass. There the largest discrepancy 
between PHOTOS EXP and KKMC was found. Bottom 
plot: The same difference in case of single photon 
runs of PHOTOS and KORALZ and ^^/i^ invariant mass 
can be seen. One can see that the source of residual dif- 
ference is in approximation of the matrix-element used 
in PHOTOS and shows the possible further improvement 
in the code. 



